■ ORIGINAL ARTICLE Print ISSN 1 738-6586 / On-line ISSN 2005-501 3 

J Clin Neurol 2014;10(3):229-235 http://dx.doi.Org/10.3988/jcn.2014.10.3.229 

CrossMark 

^HHJ <- click for updates 



Open Access 



Abnormal Brain Activity Changes in Patients with IVIigraine: 
A Short-Term Longitudinal Study 

Ling Zhao,^'^ Jixin Liu,^ Xuemei Yan,^ Wanghuan Dun,'' Jing Yang,"" Liyu Huang,^ 
Yuan Kai,^ Dahua Yu,"^ Wei Qin,^ Tian Jie,^ Fanrong Liang^ 

^School of Life Science and Technology, Xidian University, Xi'an, PR China 

"The 3rd Teaching Hospital, Chengdu University of Traditional Chinese Medicine, Chengdu, PR China 
'^Department of Medical Imaging, First Affiliated Hospital, Xi'an Jiaotong University College of Medicine, Xi'an, PR China 
^Information Processing Laboratory, School of Information Engineering, Inner Mongolia University of Science and Technology, Baotou, 
Inner Mongolia, PR China 



Background and Purpose Whether or not migraine can cause cumulative brain alterations 
due to frequent migraine-related nociceptive input in patients is largely unclear. The aim of this 
study was to characterize longitudinal changes in brain activity between repeated observations 
within a short time interval in a group of female migraine patients, using resting- state function- 
al magnetic resonance imaging. 

Methods Nineteen patients and 20 healthy controls (HC) participated in the study. Regional 
homogeneity (ReHo) and functional interregional connectivity were assessed to determine the 
focal and global features of brain dysfunction in migraine. The relationship between changes in 
headache parameters and longitudinal brain alterations were also investigated. 

Results All patients reported that their headache activity increased over time. Abnormal 
ReHo changes in the patient group relative to the HC were found in the putamen, orbitofrontal 
cortex, secondary somatosensory cortex, brainstem, and thalamus. Moreover, these brain re- 
gions exhibited longitudinal ReHo changes at the 6-week follow-up examination. These head- 
ache activity changes were accompanied by disproportionately dysfunctional connectivity in 
the putamen in the migraine patients, as revealed by functional connectivity analysis, suggest- 
ing that the putamen plays an important role in integrating diverse information among other 
migraine-related brain regions. 

Conclusions The results obtained in this study suggest that progressive brain aberrations in mi- 
graine progress as a result of increased headache attacks. J Clin Neurol 2014;10(3):229-235 
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Introduction 

Migraine is an idiopathic headache disorder defined as an epi- 
sodic disturbance primarily manifesting as head pain and touch 
sensitivity, accompanied by nausea and light sensitivity.^ Neu- 
roimaging studies of migraineurs have described functional 
abnormalities in brain regions associated with pain processing 

© This is an Open Access article distributed under the terms of the Cre- 
ative Commons Attribution Non-Commercial License (http://creative- 
commons.org/lic enses/by-nc/3.0) which permits unrestricted non-com- 
mercial use, distribution, and reproduction in any medium, provided the ori- 
ginal work is properly cited. 



as a result of repeated headache attacks, including enhanced 
cortical excitability^ and altered pain modulatory systems.^"^ 
However, the exact process underlying these brain abnormali- 
ties has not been fully elucidated, and whether or not migraine 
causes cumulative dysfunctional activities when the headache 
activity increases is still largely unknown. 

Due to frequent migraine-related nociceptive input, several 
brain regions in migraine patients (PM) have been reported to 
exhibit abnormal functional activity at rest,^ leading to irregu- 
lar brain circuits associated with pain-related information pro- 
cessing. Moreover, some studies patients have classified into 
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two subgroups according to their migraine attack frequency 
and attack history.^'^ Comparisons between the milder and 
more severe patient groups were aimed at exploring the brain 
functional changes underlying the pathophysiology of mi- 
graine progression. While these studies have to some degree 
quantitatively identified the predilection site of possible brain 
functional changes in migraine, longitudinal studies are still 
needed. Accordingly, the aim of the present study was to char- 
acterize the longitudinal brain functional abnormalities occur- 
ring in migraine within relatively short periods of time. We 
hypothesized that the presence of an increased number of mi- 
graine attacks occurred between the repeated observations 
would be indicative of longitudinal migraine-related brain 
dysfunction. 

We have shown previously that migraine has a greater influ- 
ence in females and leads to more dysfunctional brain activity 
in their resting functional networks compared to males.^ To test 
our hypothesis, we focused on a group of female migraineurs 
without aura. Regional similarity and functional connectivity 
analysis were employed to investigate both the local site of the 
brain abnormalities and the entire dysfunctional network in mi- 
graine. 

Methods 

All research procedures were approved by the West China 
Hospital Subcommittee on Human Studies and were conduct- 
ed in accordance with the Declaration of Helsinki. All subjects 
gave written, informed consent to participate after the experi- 
mental procedures had been fully explained and they were in- 
formed that they could stop participating at any time. 

Participants 

Nineteen right-handed PM without aura [all female, age 21 .8± 
2.3 years (mean±SD), migraine duration of 9.1 ±2.6 years] 
who did not have any clinical affective disorder were recruit- 
ed. Twenty, education- and gender-matched, healthy, right- 
handed healthy controls (HC; age 22.4±3.1 years) were recruit- 
ed from the local community. The controls had experienced 
no headache during the previous year and had no family mem- 
bers who suffered regularly from a migraine or other head- 
aches. All of the patients were screened in accordance with the 
International Headache Society criteria.^ The diagnostic crite- 
ria of the International Headache Society for migraine without 
aura include the occurrence of at least five headache attacks 
that fulfill the following criteria: 1) a unilateral and/or pulsat- 
ing headache, 2) headache attacks lasting 4-72 hours (untreat- 
ed or unsuccessfully treated), 3) presence of nausea and/or 
vomiting, photophobia, and phonophobia during the headache, 
4) and the headache being disabling.^ The exclusion criteria 



were 1) macroscopic brain T2-visible lesions on magnetic res- 
onance imaging (MRI) scans, 2) existence of a neurological 
disease, 3) pregnancy or menstrual period, 4) alcohol, nico- 
tine, or drug abuse, or 5) claustrophobia. 

All subjects submitted to a resting- state functional MRI 
(flMRI) scan at study entrance (i.e., baseline), and all of the pa- 
tients underwent a second MRI scan within 6 weeks (39-45 
days) using the same MR imager and acquisition protocol. For 
all of the patients, the scans were not performed within 72 
hours before, during, or 24 hours after a migraine attack. The 
drugs used for the prophylaxis of migraine were stopped 6 
weeks before the experiment; however, the subjects were al- 
lowed to take pirprofen when their headache was difficult to 
endure. Detailed information about the patients' drug intake 
was obtained. The Medication Quantification Scale was used 
to quantify the consumption of the drug at each visit, and the 
dosage and duration of the drug used are presented as scalar 
values. Patients were required to keep a headache diary to re- 
cord their headache activity, including the migraine attack fre- 
quency, migraine attack duration, and average pain intensity of 
the headache (on a scale from 0 to 10 scale, with 10 being the 
most intense pain imaginable). 

fMRI data acquisition 

All subjects underwent resting-state fMRI scans using a 3-T 
MR system (GE EXCITE, Milwaukee, WI, USA) with an 
eight-channel phased-array head coil. Routine T2-weighted 
images [repetition time (TR)/echo time (TE)=4680.0 ms/105.2 
ms, matrix size=256x256, field of view (FOV)=256x256 
mm^, and section thickness=5 mm] and Tl -weighted images 
(TR/TE=2360.4 ms/21.9 ms, matrix size=256x256, FOV= 
256x256 mm^, and section thickness=5 mm) were examined 
by an expert radiologist to exclude the possibility of clinically 
silent lesions. The fMRI images were obtained using an echo- 
planar imaging (EPI) sequence with the following parameters: 
30 continuous slices with slice thickness=5 mm, TR=2000 ms, 
TE=30 ms, FOV=240x240 mm^ and data matrix=64x64. A 
total of 205 volumes was acquired for each subject, requiring 
a total scan time of 410 s. Subjects were instructed to rest with 
their eyes closed, not to think about anything in particular, and 
not to fall asleep. 

fMRI data preprocessing 

The first five volumes were discarded to eliminate the non- 
equilibrium effects of magnetization and to allow subjects to 
acclimatize to the scanning environment. Image preprocessing 
was carried out using Statistical Parametric Mapping 5 soft- 
ware (http://www.fil.ion.ucl.ac.uk/spm). All data sets were 
processed using the following steps: 1) compensation of sys- 
tematic slice-dependent time shifts, 2) elimination of system- 
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atic odd-even slice intensity differences due to the acquisi- 
tion of interleaved images, 3) rigid-body correction for 
geometrical displacements caused by head movement, 4) spa- 
tial normalization with the Montreal Neurological Institute 
(MNI) EPI template image using an optimum 12-parameter 
affine and nonlinear cosine basis function transformation, 
resampled to 3-mm isotropic voxels, and 5) physiological 
and high-frequency noise correction by applying a band-pass 
filter from 0.01 to 0.1 Hz. No subjects exhibited head mo- 
tions exceeding 1 mm of translation or 1 degree of rotation in 
any direction. There were no significant between-group dif- 
ferences in head motion (maximal translational and maximal 
rotational head deviation). 

Regional homogeneity method 

The regional homogeneity (ReHo) method was used to mea- 
sure the intracluster information similarity, which might reflect 
the coherence of spontaneous neuronal activity. Recent stud- 
ies found that ReHo can be a highly reliable (test-retest) mea- 
sure of brain activity, and indicated that the method may 
serve as a valuable imaging feature for exploring pathological 
changes in brain function. By computing Kendall's coeffi- 
cient of concordance^^ for each voxel, the ReHo characterized 
the similarity of the time series of a given voxel to those of its 
nearest neighbors on a voxel-by- voxel basis. Individual ReHo 
maps were generated, with a Gaussian kemel with a full-width 
at half-maximum of 4 mm being used to smooth the images in 
order to reduce noise and residual differences.^ All of the ReHo 
analytical procedures were performed using the Resting- State 
fMRI Data Analysis Toolkit (X.-W. Song et al., Beijing Nor- 
mal University, Beijing, China; http://www.restfmri.net). A 
two- sample ^test was used to compare the ReHo results be- 
tween the PM and HC groups, and a paired ^test was used to 
compare changes in the ReHo values between the two resting 
scans in PM. The false discovery rate (FDR) was used to cor- 
rect the multiple comparisons. 



Functional connectivity analysis 

A network connectivity analysis was employed to assess the 
network properties in any abnormal brain regions. Building 
upon the patients' ReHo maps in the data sets, we investigated 
the distribution of the voxels that exhibited significant chang- 
es in ReHo values. The MNI coordinates of the maximally ab- 
normal voxel within an anatomical area (minimum cluster size 
of 3 voxels) were chosen, and a sphere with a diameter of 12 
mm was drawn around these coordinates. The voxels located 
in the white matter and ventricles within each sphere were re- 
moved to ensure the integrity of their structure and function. 
This process yielded several regions of interest (ROIs), which 
were considered as a set of nodes in the subsequent connectiv- 
ity analysis. The mean time courses from the deep white mat- 
ter, ventricles, and the six rigid-body motion parameters were 
regressed out from the fMRI time series of the selected ROIs. 
We computed the mean time series of each ROI and obtained 
the Pearson correlation coefficient matrix between all possible 
connections of the region pair within each individual patient, 
and all correlation values were converted using a Fisher Z- 
transformation. Nonparametric permutation tests were used to 
detect the presence of any significant group differences in the 
present study. 

Results 

Headache parameters 

The headache attack frequency was significantly higher at the 
6-week follow-up examination (9.88±5.95) than at study en- 
trance (4.76±3.27,/?<0.05) (Table 1), whereas the pain inten- 
sity (5.05±1.1 and 5.6±0.9, respectively, /?>0. 05) (Table 1) 
and headache attack duration (5.09±3.18 and 7.95±6.6 years, 
respectively, /?>0. 05) (Table 1) did not increase over time. 

Alteration in brain functional activity 

Between-group differences (i.e., PM vs. HC) in brain local ac- 
tivity were found for the thalamus (THA), putamen (PUT), 
brainstem (BS), cingulate cortex, inferior parietal gyrus [Brod- 



Table 1. Demographic characteristics of subjects 



Patients with migraine (n=19) 

Information p-value 

Time point 1 Time point 2 

Age (years) 21.8±2.3 - N/A 

Education (years) 14.7±2.0 - N/A 

Disease duration (years) 9.1 ±2.6 - N/A 

Migraine attacl<s during past 6 weel<s 

Attack duration (hours) 5.09±3.18 7.95±6.60 0.2 

Attack frequency (times) 4.76±3.27 9.88±5.95 <0.05 

Average pain intensity 5.0± 1 . 1 5.6±0.9 0.3 

Values are expressed as mean±SD. 
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Fig. 1. Between-group differences in lo- 
cal brain activity in patients and healthy 
controls. 



Table 2. Foci with significant differences in regional homogeneity 
values from patients versus healthy controls (p<0.05 FDR cor- 
rected). 



Regions 




X 


Talairach 

y 


z 


f value 


SFG 


L 


-27 


65 


11 


-4.72 




R 


36 


26 


51 


4.08 


OFC 


L 


-3 


43 


-20 


4.3 




R 


3 


46 


-20 


4.6 


ACC 


L 


-3 


17 


-8 


-4.3 




R 


3 


14 


-6 


-4.63 


MCC 


L 


-6 


-10 


42 


-4.15 




R 


6 


-10 


31 


-4.44 


Precuneus 


L 


-18 


-74 


40 


-4.47 




R 


27 


-71 


34 


-4.92 


Putamen 


L 


-12 


5 


-10 


-3.71 




R 


15 


11 


-8 


-3.91 


Brainstem 


L 


-9 


-33 


-21 


-4.07 


Thalamus 


R 
L 


-9 


-17 


1 


-4.67 


Sll 


R 
L 


-53 


-30 


35 


4.56 




R 


59 


-30 


46 


4.85 


STG 


L 


-53 


-43 


21 


4.59 




R 


42 


-55 


14 


4.37 


Occipital cortex 


L 


-21 


-64 


-7 


-4.1 




R 


45 


-87 


-1 


4.79 


Cerebellum 


L 


-21 


-42 


-26 


-5.23 




R 


27 


-63 


-37 


4.4 


ACC: anterior cingulate cortex, FDR: false discovery rate, L: left, 
MCC: middle cingulate cortex, OFC: orbitofrontal cortex, R: 
right, SFG: superior frontal gyrus, Sll: secondary somatosensory 
cortex, STG: superior temporal gyrus. 



mann's area 40], hippocampus, orbitofrontal cortex (OFC), 
and the occipital cortex ^<0.05, FDR-corrected) (Fig. 1, Table 
2). Moreover, the ReHo values for PM were significantly de- 
creased in the PUT (MNI coordinate center 27, 6, 0), BS (-3, 
-27, -3), THA (12, -15, 12), temporal cortex, and cerebellum, 
and were significantly increased in the OFC (3, 33, -15) and 
secondary somatosensory cortex (SII, BA40; -54, -45, 39) rel- 
ative to the baseline level ^<0.05, FDR-corrected) (Fig. 2). 
Significant correlations were found between changes in head- 
ache attack frequency and ReHo values in the OFC and SII 
(p<0.05X but not in the PUT, THA, and BS (Fig. 2). 

Abnormality in brain functional connectivity 

To assess the fianctional connectivity among those local abnor- 
malities, the PUT, BS, THA, OFC, and SII were chosen as the 
ROIs for the functional connectivity analysis (see Methods). 
As can be seen from Fig. 3, the changes in the functional con- 
nectivity intensity in the eight connections were significantly 
correlated with the changes in headache attack frequency 
^<0.05) (Fig. 3). However, not all of the connections exhibit- 
ed significant between-group differences. Five significantly 
increased connections were found on different resting scans 
(nonparametric permutation tests, p<0.01; red and blue lines 
in Fig. 3B). Analysis of the changes in the intensities of the 
PUT-THA and PUT-OFC connections revealed a positive cor- 
relation with the changes in headache attack frequency (red 
lines in Fig. 3B). The PUT- SII and PUT-BS connections were 
significantly negatively correlated with the changes in head- 
ache attack frequency (blue lines in Fig. 3B). 
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Fig. 2. Longitudinal brain changes in re- 
gional homogeneity (ReHo) occurred in 
migraine patients (PM) within a relative- 
ly short period of time (p<0.01, correct- 
ed for the false discovery rate). Warm 
colors indicate ReHo increases in PM 
relative to the baseline level. BS: brain- 
stem, OFC: orbitofrontal cortex, PUT: pu- 
tamen, Sll: secondary somatosensory 
cortex, THA: thalamus. 
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Fig. 3. Abnormal functional connectivity 
and its relationship with changes in head- 
ache attack frequency. The putamen 
(PUT), orbitofrontal cortex (OFC), sec- 
ondary somatosensory cortex (Sll), brain- 
stem (BS), and thalamus (THA) were 
chosen as the regions of interest for the 
functional connectivity analysis (red balls). 
A: The gray lines indicate a significant 
correlation between the changes in at- 
tack frequency and the ratio of changes 
in functional connectivity intensity (p< 
0.05). B: The color connections indicate 
the significant changes in functional con- 
nectivity intensity in the migraine group 
between the repeated observations, and 
represent significant correlations be- 
tween the changes in attack frequency 
(red lines: positive correlations; blue 
lines: negative correlations). 



Discussion 



The main results of this study are as follows: 1) migraine dis- 
rupted the brain activity in several pain-processing brain re- 
gions, 2) the changes in ReHo values in these regions between 
repeated observations were significantly correlated with the 
changes in headache attack frequency, and 3) longitudinal al- 
terations in connectivity characteristics converged at the PUT. 

Abnormal ReHo value changes were found at the follow-up 
examination in the PUT, THA, BS, OFC, and SII. Several 
studies have found a network of brain regions (pain matrix) 
that are consistently activated in response to acute painful 
stimuli, and some of these brain areas also exhibit dysfunc- 
tional activities in chronic pain conditions. ^'^^ Our findings are 
consistent with those of previous studies. The prefrontal cor- 
tex, THA, basal ganglia, parts of the BS, and the somatosen- 



sory cortex are the main regions that process and regulate pain 
information. ^^'^^ ReHo measures the similarity of low-fi'equen- 
cy fluctuations within a local brain area, which may reflect the 
coherence of neuronal activity. ^^'^^ Hence, abnormal ReHo 
values indicate that resting neural function in the aforemen- 
tioned brain regions was abnormally synchronized in PM be- 
tween the two different observation points, which may reflect 
differences in the pain-related information processing induced 
by the increased headache activity. The significant correlation 
between the changes in headache attack frequency and ReHo 
values may support our assumption. Several studies have 
shown that chronic headache attack can be understood not only 
as a dysfunctional perceptual state but also as a result of brain 
plasticity in response to long-term chronic pain.^^ Hence, these 
brain abnormalities may be the result of aberrant develop- 
ment of the brain in the presence of ongoing central changes. 
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rather than as an adventitious insult or lesion to a normal 
brain. This process may be maladaptive and directly related to 
the worsening of headache activity in migraine. The current 
results suggest a relationship between progressive brain func- 
tional alteration and increased migraine attack frequency in 
patients. However, whether these longitudinal ReHo changes 
can be employed to predict progressive neuronal changes or to 
evaluate the worsening of clinical symptoms in particular pa- 
tients requires further investigation. 

The BS is the main station for integrating the cognitive and 
autonomic processing of nociceptive information,^^ and is con- 
sidered to be made dysfunctional by migraine pathophysiolo- 
gy 20,21 Positron-emission tomography has highlighted neuro- 
nal activation in the BS in PM.^^'^^ Similar neuronal changes 
were found in experimental animals during trigemino vascular 
activation.^^'^^ Moreover, measurement of the neural activity 
before and after successful treatment in migraine revealed that 
the BS remained activated,^^'^^ indicating that the BS might 
have generated the migraine attack. ^^'^^ The findings of the 
present study cannot be used to determine the cause-and-effect 
relationship between BS modulation and migraine symptoms. 
The longitudinal brain changes in the BS were independent of 
changes in the attack frequency. We assumed that changes or 
dysfunction in the BS may present a migraine-related brain 
state, creating a referred headache, which caused the corre- 
sponding dysfunction of the brain regions involved in the pain 
processing networks. 

Another important finding of this study was the identifica- 
tion of the disproportionately dysfunctional connectivity in the 
PUT, while acting as a critical intermediate station for infor- 
mation processing among the selected ROIs in migraine. Sev- 
eral brain imaging studies have found abnormal activation in 
the basal ganglia, and notably in the PUT,^^ in migraine. More- 
over, the gray matter in the basal ganglia was altered in mi- 
graineurs with a higher attack frequency and longer duration 
of disease.^ One recent study found significantly lower activ- 
ity in the PUT for the response to painful heat stimulation in 
high frequency episodic migraine as compared to low fi-equen- 
cy episodic migraine.^ That finding may support the inference 
that the dysfunction in the PUT would be aggravated by an in- 
creasing migraine frequency. Our results are consistent with 
those obtained from previous studies. We found that the PUT 
serves to integrate diverse information among other migraine- 
related brain regions, and we infer from this that preventative 
treatments that are designed to target the dysfunction of the 
PUT in migraine may slow the disease process and relieve 
pain in migraineurs. 

Two issues raised by the present study need to be addressed 
further. First, little is known about the differences between 
progressive patients and nonprogressive patients with respect 



to dysfiinctional brain organization in migraine pathophysiol- 
ogy. Second, the interregional functional connectivity analysis 
was not corrected for in our multiple comparisons. The com- 
plex patterns of negative and positive correlations between 
functional connectivity and migraine attack frequency were 
difficult to explain; these results should thus be interpreted with 
caution. 

In summary, a sustained increase in headache attack fre- 
quency and longitudinal functional brain activity changes 
were found to occur in migraineurs within a relatively short 
period of time. These findings may reflect brain alteration dy- 
namics resulting from a developed headache. If these results 
can be replicated in a large sample, it may help us to better un- 
derstand the progression of migraine and improve our ability to 
implement effective treatments. 
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